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ABSTRACT 


The one thing among many that must be prepared for every flight is the making 
of a flight plan. Pilots must use charts or graphs from the flight manual to compute the 
fuel flow that is essential to a flight plan. Since this requires many steps of 
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manual, it is time consuming and increases the probability of making a mistake. This 
problem cobstructs the execution of various mission changes and continuous sorties. 

A computer program for personal computer or hand-held calculator is developed 
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I. INTRODUCTION 


A. BACKGROUND 

A pilot must be furnished with flight performance information for his aircraft 
that will permit him to operate the aircraft efficiently at different airspeeds, different 
altitudes, different values of gross weight and different external loadings. 

This information is collated from a combination of theoretical determinations 
and flight test data and presented to the pilot in a flight handbook in either a tabular 


format or in a series of graphs, or in combinations of these methods. 


B. FUEL FLOW RATE 

Probably the most important parameter in optimal performance of an aircraft is 
the rate at which fuel is consumed. It 1s obvious that when the fuel is exhausted the 
aircraft will no longer fly, so the pilot is very interested in the rate at which the fuel is 
being consumed. But also of importance is the fact that as the fuel is burned the weight 
of the aircraft is decreased and the performance of the aircraft is changed. It 1s 
therefore quite important that the pilot be able to forecast the rate of fuel consumption 
for each portion of each flight. 

Fuel flow rate for a turbo-jet aircraft is a direct function of the thrust required 
for the specified flight condition, and in equilibrium flight thrust required is equal to 
the total drag of the aircraft. 


C. TOTAL DRAG 

The total drag of an aircraft in equilibrium flight may be subdivided into induced 
drag (the drag associated with the generation of lift) and parasite drag (the drag 
associated with the resistance of the fuselage to motion through the air). The induced 
drag is a function of the altitude of the aircraft (air density) and the weight of the 
aircraft, and the parasite drag is also a function of the altitude as well as being a 
function of the profile drag of the aircraft. The aircraft drag, which varies with the 
configuration and external loading (weapons, fuel tanks, et cetera) 1s usually expressed 
in terms of a Drag Index (DI), with the zero Drag Index usually being a measure of 


the drag of the aircraft in the ‘clean’ configuration. 
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D. THRUST 
Although it is the thrust of the engines that provides the performance for a 
turbo-jet aircraft, as long as there is sufficient thrust available to provide the required 


performance, the pilot is more interested for the flight planning in the fuel flow rate. 


E. FLIGHT PLANNING 

To plan a flight, the pilot must determine what will be the loading of the aircraft 
(gross weight) and what effect this loading will have on the drag (Drag Index). The 
pilot then enters the flight handbook to determine the fuel flow rate for the airspeed 
and altitude at which he desires to fly. 

If the gross weight, Drag Index and airspeed are exactly the same as values listed 
in the handbook tables or exactly the same as the valued on the graphical lines, the 
determination of the fuel flow rate is simply a matter of reading (correctly) the values 
from the information present in the handbook. 

If, however, any of the input parameters are not identical to the handbook data 
points, the pilot must interpolate between the values presented. If weight, Drag Index, 
airspeed and altitude are all off-presentation points, the pilot may be forced to do 
several cross interpolations. This may even involve interpolations between values listed 
on different pages. Not only is this a time consuming and difficult task, but the 


possibility for error is expanded with each interpolation. 


If. APPROACH TO THE PROBLEM 


A. STATEMENT OF THE PROBLEM 

A method was sought that would permit the pilot to rapidly obtain the fuel flow 
rate for various flight combinations of gross weight, Drag Index and airspeed. And, 
regardless of the ease of access to the information and completeness of the data 
obtained, this method must provide fuel flow rate at the same degree of accuracy as the 
flight handbook. 

In addition, it was planned that the method of data retrieval should be as simple 


as possible with an ultimate goal of using a hand-held computer for the solutions. 


B. BASE LINE APPROACH 

The standard method of using the flight handbook is for the pilot to search and 
follow several lines and tables of the handbook to determine the required fuel flow rate. 
This requires that altitude, gross weight and Drag Index must be considered, but this 
method 1s too long, takes too much effort and 1s open to errors. 

The procedure can be simplified by computerizing the flight manual. Most 
attempts at this method have involved the development of curve-fitting routines that 
will-supply numerical modeling of the handbook data. Although the end use of these 
equations is quite rapid, the development of these routines is a long and arduous task. 
Inasmuch as each equation usually fits only a small part of the data, the computer 
programs for such a method are usually quite lengthy. 

The approach undertaken in this project was the use of the equations for an 
‘IDEAL’ aircraft to develop simplified equations that can be used in a short computer 
program. 

Three steps were undertaken in the development of these equations: 

a) The fuel flow rate data in the handbook is converted to thrust required. 


b) Ihe perust required equations are adjusted so that they represent the actual 
aircraft. 


c) The corrected thrust required equations are then reconverted to flow rate as the 
program output. 
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C. DETAILED METHODS 
1. Concept of K, and K, 
Determinations are made for equilibrium, unaccelerated conditions where drag 
equals thrust. 


The basic drag equations for an ideal aircraft are used [Ref. I:p. 160] 


Ch — De aa Cp, 
Cp, = C)7/me AR 
D = (1/2)p V7 SCp 


All of the principal items of flight performance involve steady state flight 
conditions and equilibrium of the airplane. For the airplane to remain in steady level 
flight, equilibrium must be obtained by a lift equal to the airplane weight. And the 
airplane drag equals the thrust required to maintain steady level flight. 

The total drag of the airplane is the sum of the parasite and induced drag. 
Parasite drag is the sum of pressure and friction drag which is due to the basic 
configuration and, as defined, is independent of lift. The equation of parasite drag can 


be expressed as a function of squared velocity. 
Dy = (1/2) pS v? 
Let K, = Sy pS 
Dee Ke 


Induced drag is the undesirable but unavoidable consequence of the 
development of lift. 

In the airfoil, the local lift vector is aligned perpendicular the the local relative 
wind and hence is inclined behind the vertical by the angle @ . Consequently, there is a 
component of the local lift vector in the direction of Vg ; that 1s, there is a drag 
created by the presence of downwash. 

This drag is defined as induced drag, denoted by D; [Ref. I:pp. 152-155] 


Z 2 
Dj = (1/2)p S V* Cp, 
Cp, = C7 / (me AR) 


D; = (C,)* /(2me AR)p S V" 
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= [(2W/p S V*)* / (2m e AR) p S V2 
= (4W? pS V*)/ (p? S* V4 we AR) 
= ((2W2) / (pS we AR)I(1/V2) 


If the whole item of mght side except the squared velocity is changed to some 


constant, this equation can be expressed function of squared velocity. 
Let Ky = (2W*)/ (pS me AR) 
D. = K4(1/V7) 


While the parasite drag predominates at high speed, induced drag 
predominates at low speed. Figure 2.1 illustrates the variation with speed of the 
induced, parasite, and total drag for a specific airplane configuration in steady level 
flight. 

The curve of drag or thrust required versus velocity shows the variation of 
induced, parasite, and total drag. Induced drag predominates at low speeds when the 
airplane is operated at maximum lift-drag ratio, (L/D),,,, the total drag is at a 
minimum and the induced and parasite drags are equal. 

The effect of a change in weight on the thrust required is illustrated by Figure 
2.2. The changes in the drag curves with variations of airplane weight, configuration, 
and altitude furnish insight for the variation of range, endurance, climb performance, 
etc., With these same items. 

First, the primary effect of a weight change is a change in the induced drag at 
any given speed. Thus, the greatest changes in the curves of thrust and power required 
will take place in the range of low speed flight where the induced effects predominate. 
The changes in thrust and power required in the range of high speed flight are 
relatively slight because parasite effects predominate at high speed. The induced effects 
at high speed are relatively small and changes in these items produce a small effect on 
the total thrust required. [Ref. 2:p. 99] 

Second, Figure 2.3 illustrates the effect on the curve of thrust required of a 
change in the equivalent parasite area f of the configuration. Since parasite drag 
predominates in the region of high flight speed, a change in f will produce the greatest 


change in thrust and power required at high speed. Since parasite drag 1s relatively 
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small in the region of low speed flight, a change in f will produce relatively small 
changes in thrust required at low speeds. The principal effect of a change in equivalent 
parasite area of the configuration is to change the parasite drag at any given airspeed. 
US 23a, TO 

Third, a change in altitude can produce significant changes in the curves of 
thrust required. The effect of altitude on these curves provide a great part of the 
explanation of the effect of altitude on range and endurance. [Ref. 2:p. 101] 

Figure 2.4 illustrates the effect of a change in altitude in the curves of thurst 
required for a specific airplane configuration and gross weight. As long as 
compressibility effects are negligible, the principal effect of increased altitude on the 
curve of thrust required is that specific aerodynamic conditions occur at higher true 
airspeeds. 

As being investigated above, the total drag, same as thrust required, is the 
function of squared velocity with different constant. 


Total drag equation is 
D = K, V* + K/ Vv? 
Where Ky, = Spy. pS 
Ky = (2 W*)/(p S me AR) 


The first step of this project is to determine the K, and K> values to obtain 
the thrust required. It can be seen that the K, is the function of altitude and drag 
index, K4 is the function of gross weight and altitude. So if K, and K> values can be 
computed with any flight condition, the thrust required will be solved at any airspeed. 

2. Determining the K, 

To determine the K), use the minimum drag concept. 

At the minimum drag point, the airplane can be flown maximum endurance 
and the induced drag and parasite drags are the same. 

The first step for determine the K, is to investigate the thrust equation. Ihe 
thrust equation can be established with fuel flow for each airspeed and some factor C. 

Pawel = PURE PLOW x C 
Since the fuel flow of the turbojet powered airplane is proportional to the 


airspeed, the thrust required is proportional to the fuel flow. 
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The factor C 1s the value of the military power devided by fuel flow at military 
thrust. 

C = MILITARY THRUST / FUEL FLOW ( at military thrust ) 

The military thrust and fuel flow information for any condition can be 
obtained from NATOPS FLIGHT MANUAL. With this factor C, the thrust required 
for each speed can be computed by fuel flow,came from NATOPS FLIGHT 
MANUAL, multiplies factor C. Table 1 and Figures 2.5 and 2.6 show the relationship 


and the result between thrust and velocity. 


TABLE 1 
THRUST REQUIRED AND FUEL FLOW VERSUS VELOCITY 


VELOCITY(KTAS) FUEL FLOW(LBS/H) THRUST REQUIRED(LBS) 
360 7,725 6,558 
400 9014 7652 
440 10°623 9'018 
480 12°462 10°579 
520 14°652 127438 
560 Sa 14°752 
600 217393 187160 
MIL 25°680 21800 


Example) military thrust  ; 21800 lb 
fuel flow (at mil) ; 25680 Ibs/hr 
Factor C = 21800/25680 
= 0.8489] 
Thus 
Thrust required = 0.84891 x Fuel Flow 


Maximum endurance speed at specific gross weight, altitude, and drag index 
can be figured out with the graph of maximum endurance mach number in NATOPS 
MANUAL Figure 2.7. This true mach number needs to be transfered to KTAS to 
substitute into the drag equation. 

Next step is to compute the minimum drag. 

Minimum drag can be computed from the fuel flow at maximum endurance 
multiplied by the factor C. The maximum endurance fuel flow can be obtained from 
the NATOPS MANUAL Figure 2.8. 

And so the minimum drag is the twice of the parasite drag. 

D = 2Da= 2D aa 


total — min 
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EFFECTIVE GROSS WEIGHT - 1000 POUNDS 


TOTAL FUEL FLOW — 1000 POUNDS PER OUR 


Thus the drag equation can be expressed. 
= = 2 
|) = Be = O52) Coy 


where p = density at given altiude 
V = velocity at minimum drag 
S = 


wing area of specific airplane 


[D. can be obtained from the above equation, because all values except oF are 
known, and the drag equation is a function of squared velocity with the known 0), 
The total value of right side except V~ is the desired constant Ky. 
3. Determining the K, 

The procedure to compute K- is just similiar to that for computing Ky. ft 1s 
convenient to use the same minimum drag concept. At the minimum drag point, the 
minimum thrust required is twice the induced drag amount. 

Use the same maximum endurance velocity and minimum thrust required. 


The induced drag equation is 


D: = ee 2 


1 


= (2W’) / (we AR p S\(1/V2) 


where w = Gross weight of specific airplane 


= Wing area of specific airplane 


Air density at given altitude 


p 


v = Maximum endurance velocity 


Since the minimum thrust required and all items of the right side are known 
value except 1/7 e AR, the value of 1/m e AR can be solved easily. After obtained the 
1/tm e AR value, substitute this value into the above induced drag equation, and 
calculated the numerical value of the right side except the 1 iv. This numerical value is 
just the desired constant K value. 

But if the raw K, and K>4 values are used to compute the thrust required, 
there will exist a deviation of more than 10 percent at high speed. Therefore one needs 
to modify the K, and K, to fit the curve to the FLIGHT MANUAL data. 

The fitting method used in this project is the trial and error method. The 
Table 2 and Figure 2.9 illustrates the difference between the modified and raw value of 
Ky, Ko. 
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4. Relationship between K,, K, and Variation Factors 

Next consideration is to make the equations of K, and K, functions of drag 
index, altitude and gross weight. 

Since the Ky is the function of drag index and air density, K, will be changed 
by different drag index and different altitude. Thus it is needed to develop a 
relationship of K, with the drag index and altitude. 

To accomplish this, it is convenient to select some reference altitude and gross 
weight, because the relationship of K, and drag index is independent of altitude and 
gross weight. Each K, value for every drag index can be calculated with the previous 
procedure. After calculation of all of the Ky values, a graph of K, versus drag index is 
plotted to find the relationship. 

In this project, the altitude and gross weight are selected to standard sea level 
and 40,000 lbs. The final values and relationships are shown in Table 3 and Figure 
2.10. 


WABLES 
DRAG INDEX VS K1 


Drag Index kK, 

0 0.011534 

20 OO 3835 

40 0.014807 

60 0.016780 

80 0.018480 

100 0.020527 

120 0.022361 


Fortunately, as can be seen in the result of the graph, the relationship between 


K, and the drag index is linear. So the equation K, as a function of drag index is 
K, = (0.0115 + 0.0001 x Drag Index) 


In this equation, the most important thing is the slope of the graph, because 
the starting value of K, was modified to fit the published data. 
The K, equation with modified starting value of Ky 1s 


K, = (0.014 + 0.0001 x Drag Index) 


This equation was verified by calculating the drag at each speed. The results 


are shown in Tables 4, 5, 6. 
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TABLE 4 
THRUST REQUIRED DEVIATION lee: 


) ALTITUDE GROSS WEIGHT | 40,000 LBS Fp.r | 20 


THRUST 21,800 LB 
MIC Taher THRUS® 
FUE PSE LOW -2o,680 LB/H 


AIRSPEED] FUEESELOW | PARUST (LB)| THRUST sees 







FUEL FLOW x 0.8489 









DEVIATION 












CALCULATION 





(ers) CBS i) HAND BOOK 


Ie See 


FUEL FLOW 5,600 (Leis ey 


MAXIMUM ENDURANCE THRUST (min) 4,754 Cols) 


ATR SPEED (KTAS ) 


E> 0.014+0.0001x20 = 0.016 
VARIATION OF K ae a: 


From Tables 4, 5, 6, the following can be observed. 
First, the thrust required, i.e., calculated drag, is almost same as published 


data at low drag index range, but the thrust required is much less than published data 








at high drag index range. 
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eee 
THRUST REQUIRED DEVIATION AT DI 60 


pee TCDS GROSS WEIGHT | 40,000LBS bp.r] 60 
THRUST 21,800 LB 

MILITARY Tromeot SH hUEL FLOW x v.59 25 
FUEL FLOW | 24,425 LB/H 


eee ee) FUEL PLOW |@PHRUST (LB)}| THRUST (LB) DEVIATION 














PEBS/ rH) HAND BOOK CALCULATION 


PUEL FLOn Gy 200 (LBS/H) 
MAXIMUM ENDURANCE TARUSD amin) a, 504 (LB) 
ine SPEED 240 (KTAS) 
0.014+0.0001x60 = 0.02 
VARIATION OF K 2 ; 


Second, in all drag index ranges, the thrust required is much less than 
published data at high air speed range (approximately above Mach number 0.8). 
To compensate for the first problem, that is low thrust required in high drag 


index range, it is needed to increase the K, value, but the thrust required value is 
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TABLE 6 
THRUST REQUIRED DEVIATION AT Digs? 


EY LIMICLIO © GROSS WEIGHT |40,000 LBS ale 
THRE 21,800 LB 

MILITARY THRUST 
FUEL FLOW | 23-785 LB/H 


AIR SPEED} FUELMELOW PetRRUST “(LB ) Seis ere 





FUEL FLOW x9.9165 










DEVIATION 
(%) 





















(IRBisy i) HAND BOOK CALCULATION 





aoe | neal 
a Oe 
A ee 
ee 


FUEL FLOW 6,650 ( EBisy i) 


MAXIMUM ENDURANCE | THRUST (min) 6,003 (LB) 

AIR SPEED (KTAS) 
aan 0.014+0.001x80 = 0.022 7 

VARIATION OF K 


slightly larger than the published data. The way to compensate in both the low and 


high drag index range is to increase the slope of the K, equation and to decrease the 





Starting value. 
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The modified equation of K 1s 
Ky, = (0.0135 + 0.000132 x Drag Index) 


And the change of K, equation can be illustrated in Figure 2.11. 

The second problem, that the thrust required is much less than the FLIGHT 
MANUAL data in the speed range higher than Mach number 0.8, can be illustrated as 
the airplane reaches very high flight speeds, the drag rises in a very rapid fashion due 
to compressibility. Since the generalized equation for parasite drag does not account 
for compressibility effects, the actual drag rise is typified by the dashed line in Figure 
Dele. 

This phenomenon can be investigated from the viewpoint of drag coefficient. 
The lift and drag coefficients vary with the Mach number, however, the drag coefficient 
for a fixed value of lift coefficient does not change until the drag-divergence Mach 
number is reached, at which time the drag coefficient increases sharply. If wind-tunnel 
tests of a model of the airplane to be analyzed have been conducted, curves of drag 
coefficient against Mach number for various values of lift coefficient will be available 
as shown in Figure 2.13. 

Drag data above the divergence Mach number cannot be calculated with 
precision, and to achieve any accuracy at all, wind-tunnel tests or other empirical data 
must be relied on. It is noted that even wind-tunnel data may not be reliable in the 
transonic range, and therefore performance guarantees in the transonic range, based 
purely on wind-tunnel data and/or calculations, cannot be considered completely 
dependable. 

In most instances, the lift-drag relation may be expressed in an analytic form 


such as 
~ 2 
Ch = By, + Cr / te AR 


in which case performance characteristics may be analytically determined if similiar 
analytical expression for power available can be written. 


The equation 
2 
Gee Clee C.7/me AR 
D Dy loa 


applies to most aircraft through the Mach number range of 0 to about 4 and a Cy 


range from 0 to 0.6, for both trimmed and untrummed conditions, however, Ch ande 
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AIRPLANE DRAG VS VELOCITY 


i 
AIRPLANE DATA 


f = 4.3 SQFT. 
w= 12,000 LBS. 
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Figure 2.12 Typical Airplane Drag Curve. 


become functions of Mach number above Mach numbers of about M=0.8. For 
supersonic aircraft the value of “be increases by a factor of 2 to 3 through the 
transonic region and thereafter easel exhibits a slight decrease. [Ref. 3:pp. 275-278] 
The span efficiency factor e shows a steady decrease from M=0.8 to about M=2.0 
and thereafter remains almost constant at about 50 to 60% of its low speed value. 
From M=0.8 to M=1.2 the parasite drag curve quite closely resembles a sin? curve, 


and for approximation analysis may be represented by the following equation 


25) 


DRAG COEFFICIENT IN TRANSONIC 


O46 | | | 
| Drag coefficient vs. M | 





“COCR Ne 
“|_| TT | oe 8 ae -NO 
| See 





Figure 2.13 Variation of Cd in Transonic region. 
Ch =Cp +(Cp -Cp_) sin? [(M-0.8)/0.8 
Do Deg ( Dom De [C )/0.8] 
where Cp | = maximum value of parasite drag coefficient 
om 
Cp = minimum low speed value of parasite drag 


00 
coefficient. 
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The variation of e may normally be fitted by the polynomial 


e = e,-e)(M-08)7 + e,(M-0.8)° 
where e, = low speed value of e 


€5, €3; = constants which depend on the airplane. 


With the above theory and parasite drag equation, it is possible to make the 
drag equation between the Mach number M.8 to M 1.2. 


= D 
Dy = (1/2) PSCp V 
at above Mach number M.8, where Mach Drag exist, 
D. =(1/2)pS{Cy, +(C -Cy_ )sin? {(M-0.8)/0.8} ] V2 
5 = Te Doo ( yee Doo? (1 )/0.8} | 
Dot ieeaiKew) C € AG sin? {(M-0.8)/0.8} ] V2 
5. Se 1/ pa Byer Det! {( )/0.8} ] 
=K[1+ (Cp -Cp Cp sin? {(M-0.8)/0.8} } V2 


Since K, and K> are function of air density, it must be considered altitude 
correction factor of the changing altitude. [Ref. 4] 
First consideration is the lapse rate and vertical temperature structure for the 
standard atmosphere. 
The adopted value for the lapse rate(-dT/dH = a) in the troposphere is the 
following constant 
a) In terms of standard geopotential meters 
a = 0.0065°C m! 
b) In terms of standard geopotential feet (a) 
a = 0.00356616°F ft"! 
c) In terms of the c.g.s. unit of H 
a = 6.628155 x 10°° °C cm™ sec? 
Accordingly in the troposphere 
=a (b) 
where a,H and le are in consistent unit 
Second consideration is the relationship between pressure and_ vertical 


displacement. 


oy 


The following relationships are to be understood as expressed in any system of 
consistent units. 


The air 1s assumed to obey the perfect gas law, which, in the c.g.s system of 
units, may be written as 


p = (1/R)P/T) 


where 


(c) 


p density of air, gm cm? 

p pressure, dynes cm”? 

T absolute temperature,°K 

R gas constant for 1 gram of dry air, ergs om | (eK}"! 


And the air is assumed also to be in gydrostatic equilibrium and to satisfy the 
differential equation 


dp = -gdz (d) 
where, in c.g.s units, 
P pressure, dynes em”? 


p density (specific mass), gm cm”? 


Z, vertical distance, cm 


g gravitational acceleration, cm sec’? 


The vertical displacement is herein expressed in units of geopotential. 
Geopotential is defined in differential form by the equation 
G dH = gdz 


where 


(€) 


Z altitude measured positively upwards at a point 

g positive (absolute numerical) value of the acceleration 
due to gravity at the point 

H geopotential at the point 

G dimensional constant, the amount of which determines 
the magnitude of the unit of H in terms of length and 


time.(the dimension of G are in units of gZ per unit 
of H.) 


Substituting equation (b) in equation (c), equation (e) in equation (d), and 
combining the results gives 


dP/P = G/R (-dH/T, - aH) 
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= Cain a aH)/(T, - aH) 
Intergrating the right-hand member between limits of 0 and H gives 


In( EG) as SEV. 15 1 ee ail lie) 
Lie n = G/aR 
Then 
Pelee alle) = Gh) (f) 


From equation (f), using the consistent numerical values of G,a,and R 
n = 5.2561 (dimensionless) 
The Figure 2.14 shows the relationship between temperature, pressure, density 
and altitude. 
With several equation investigated in above, the useful relationship between 
temperature, pressure change and vertical displacement. 


For this project, the useful relationships are written as 


T/T. = 1-6875x 10° H 

Bie = (6.875 x 10° H >?! 
T = 518.688 x (1-6.875x 10° H) 
a? = 1.4x1714.87xT 


O 


Finally the relationship between density of air and vertical displacement is 


p/P, = (P/RT)(Po/RoTo) 
=(P/P, T/T) 


And it can be seen by computing several cases that the effect of changing 


altitude becomes larger than normal density changing ratio. 


Let Ko =(p/p,)! 


Then, K, is proportional to the Kg and K,> is inverse proportional to the Kg . 


As was shown in previous equations, K is proportional to the squared gross 
weight, and the effect of changing gross weight becomes larger than normal gross 
weight changing ratio, same as altitude change. Thus, if there was a specific reference 
value of K~ for specific gross weight, the relationship between arbitrary gross weight 


and reference gross weight can be made easily as 
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y2.7 


Ky = Ky (W/W 
eg ree 


With given three specific conditions,i.e. altitude, gross weight and drag index, 
the thrust required can be computed for any airspeed condition with equations shown 
above. K, and K>5 values can be computed with specific altitude, gross weight and drag 
index, and the thrust required for any air speed can be computed with K, and K>. 

The final goal of this project is to compute the fuel flow rate for some given 
condition, 1.e. specific altitude, gross weight, drag index and air speed. The fuel flow is 
thrust required devided by factor C. The factor C is the value that military power 


devided by fuel flow of that time, as was shown in previous step. 


C = MILITARY THRUST / FUEL FLOW (at military thrust) 
5. Factor C 
To compute factor C, it is necessary to find the relationship military thrust 
change for different altitude and fuel flow change. This is done in two steps. 
a. Computing the military thrust 
Altitude is one factor which strongly affects the performance of the turbojet 
engine. An increase in altitude produces a decrease in density and pressure and, if 
below the tropopause, a decrease in temperature. If a typical nonafterburning turbojet 
engine is operated at a constant RPM and true airspeed, the variation of thrust and 
specific fuel consumption with altitude can be approximated from Figure 2.14. The 
variation of density in the standard atmosphere is shown by the values of density ratio 


at various altitudes. Typical values of the density ratio at specific altitudes are as 


follows 
Altitude, ft: Density ratio 
SIGE [ENG) | one ee nn rr 1.0000 
RE EM ace sincscon ester ceccatvsccsscduedee 0.8617 
LOR OO OR. :... SORE... I... ORT) 
22 COCs s. css sv sseseeaveaens 0.4976 
eG Ue aa sarc v cele eduave cauiaaanw raters: 0.3099 
ALO (CC DiS) Sener en te 0.2462 
SIO ere Rn eee ee 0.1532 


If the fixed geometry engine is operated at a constant V (TAS) in subsonic 


flight and constant N (RPM) the inlet velocity, inlet ram, and compressor pressure 
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ratio are essentially constant with altitude. An increase in altitude then causes the 
engine air mass flow to decrease in a manner very nearly identical to the altitude 
density ratio. Of course, this decrease in mass flow will produce a significant effect on 
the output thrust of the engine. Actually, the variation of thrust with altitude is not 
quite as severe as the density variation because favorable decreases in temperature 
occur. The decrease in inlet air temperature will provide a relatively greater combustion 
gas energy and allow a greater jet velocity. The increase in jet velocity somewhat 
offsets the decrease in mass flow. Of course, an increase in altitude provides lower 
temperatures below the tropopause. Above the tropopause, no further favorable 
decrease in temperature takes place so a more rapid variation of thurst will take place. 
The approximate variation of thrust with altitude is represented by Figure 2.15 and can 


be computed with the equation as follows 


i eel 7) / g(1/2) 


avail — 


Whletes 7 Ob = alge) ues) 
© = Patt / Pssi 


Some typical values at specific altitudes are as follows 


Altitude, ft: Ratio of (Thrust at altitude/thrust at SSL) 


Sea level’ .......scccssetoess sean eae 1.0000 
5 OOO. .......sda alae erearen en lantemreetea mies 0.8960 
10,000 oo. ts ccsesssdenessoss os teeiiennenss eae 0.8003 
15,000 sceescstecccaces erence eae erate tees 0.7114 
20,000 .......ssscissescnntaseneeenete seems 0.6296 
25,000" sasicssassis ea eederecen a rmreneee aoreer 0.5544 
30,000 ss. cccssccsnaeeeeeeeeeeeeene rs cetee sea waas 0.4855 


Since the change in density with altitude is quite rapid at low altitude, turbojet 
takeoff performance will be greatly affected at high altitude. Also note that the thrust 
at 35,000 ft. is approximately 39 percent of the sea level value. 

The thrust added by the afterburner of a turbojet engine is not affected so 
greatly by altitude as the basic engine thrust. The use of afterburner may provide a 
thrust increase of 50 percent at low altitude or as much as 100 percent at high altitude. 
[Ref. 2:pp. 119-121] 
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When the inlet ram and compressor pressure ratio is fixed, the principal 
factor affecting the specific fuel consumption is the inlet air temperature. When the inlet 
air temperature is lowered, a given heat addition can provide relatively greater changes 
in pressure or volume. As a result, a given thrust output requires less fuel flow and the 
specific fuel consumption is reduced. While the effect of altitude on specific fuel 
consumption does not compare with the effect on thrust output, the variation is large 
enough to strongly influence a typical variation of specific fuel consumption with 
altitude. Generally, the specific fuel consumption decreases steadily with altitude until 
the tropopause is reached and the specific fuel consumption at this point is 
approximately 80 percent of the sea level value. 

Above the tropopause the temperature is constant and altitude slightly 
above the tropopause causes no further decrease in specific fuel consumption. Actually, 
altitudes much above the tropopause bring about a general deterioration of overall 
engine efficiency and the specific fuel consumption begins an increase with altitude. 
The extreme altitudes above the tropopause produce low combustion chamber 
pressures, low compressor Reynolds Numbers, low fuel flow, etc, which are not 
conducive to high engine efficiency. 

The available military thrust for this project with various altitude 1s shown 
onetaole 7 and Figure 2.15. 


TABEE 7 
MILITARY THRUST FOR DIFFERENT ALTITUDES 


Altitude (ft) Military thrust (lbs) 

Sea level 21,800 

000 19,533 
10,000 17,446 
15,000 15,508 
20,000 13.725 
25,000 12,086 
30,000 10,584 


The appropriate equation for the curve is as follows 
THRUST = 21,800 - 0.4 x ALTITUDE 


b. Computing the fuel flow at military thrust 
There are three effective factors required to change the fuel flow. One is 


altitude, another is drag index, and the other is gross weight. To get an appropriate 
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fuel flow equation as a function of altitude, drag index and gross weight, find the 
relationship between fuel flow and each effective factor, then combine the relationship 
totally. 

(1) Relationship between fuel flow and altitude. 

According to the NATOPS FLIGHT MANUAL data, the fuel flow is 
decreasing gradually with increasing the altitude with other conditions constant. To 
find a relationship between fuel flow and altitude, specific drag index and gross weight 
are fixed, to select the available data. Table 8 and Figure 2.16 show the relationship 
when drag index and gross weight are fixed to zero and 40,000 lbs. 


After curve-fitting, the fuel flow equation stated as follows 


ie = 26,000 - 0.64 x Altitude (at below 8,000 ft) 
Fy = 20,840 - 0.37 x (Altitude - 8,000 ) (at above 8,000 ft) 
where 


F, = Fuel flow considered with different altitude 
(2) Relationship between fuel flow and drag index. 
According to the NATOPS FLIGHT MANUAL data, the fuel flow is 
decreased as the drag index increased. As was considered in step (1), a specific altitude 


and gross weight will be fixed to find the relationship between fuel flow and drag index. 


Table 9 and Figure 2.17 show the relationship when altitude and gross 
weight are fixed to zero and 40,000 lbs. 

To find the relationship between fuel flow and drag index with different 

altitude, the computed fuel flow from step (1) must be used. With above two variable 


coefficients, the fuel flow equation is 


F. = F, - 16.65 x Drag Index (at below Drag Index 20) 
F. = (F, + 333.3) - 22.1 x Drag Index (at above Drag Index 20) 
Where 


F, = Fuel flow came from step (1) 


= Fuel flow considered with different altitude and 


as 
) 
| 


drag index 


(3) Relationship between fuel flow and gross weight. 
According to the NATOPS FLIGHT MANUAL, the fuel flow 


decreased a small amount with increasing gross weight. To find the relationship 
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TABLE 8 
PUEDE EO WerOeDIPFERENT ALTITUDES 


Altitude (ft) Fuel Flow (lbs/hr) 
Sea level 26,000 
4,000 25,520 
§ 000 20,840 
12,000 19,000 
TABLE 9 


FUEL FLOW FOR DIFFERENT DRAG INDEXS 


Drag Index Fuel Flow 

) 26,000 
20 25,680 
40 25,110 
60 24,425 
80 385 
100 28.225 
120 22,745 

TABLE 10 


FUEL FLOW FOR DIFFERENT GROSS WEIGHTS 


Grossweight (Ibs) Fuel Flow (lbs/hr) 
35,000 26,000 
40,000 26,000 
45,000 25,995 
50,000 25290 
55,000 25,985 
60,000 25,980 


between fuel flow and gross weight, a specific altitude and drag index must be selected 


as in the previous steps. Table 10 and Figure 2.18 show the relationship when altitude 


and drag index are fixed to standard sea level and zero. 


It can be seen from Figure 2.18 that the relationship between fuel flow 


and gross weight is linear. The final equation of fuel flow with different three variables, 


1.e., altitude, drag index and gross weight, can be written as follows 


= F, -0.002x(Gross weight-40,000) (at heavier than 40,000Ibs) 


(at lighter than 40,000 Ibs) 


where 


F, = Fuel Flow came from step (2) 
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ey Fuel Flow considered with different altitude, 


drag index and gross weight 


Now the military thrust and fuel flow for any altitude, drag index and 
gross weight can be computed with the above step (1) and (2). The desired factor C can 
be computed with the military thrust and fuel flow. 


Factor C = Military thrust/ Fuel Flow 


To compute the fuel flow for any various specific condition is the final 
goal of this project. The fuel flow can be computed with thrust required, 1.e. total drag, 
devided by factor C. The total drag and factor C had been computed in previous step 
already. 


FUEL FLOW = THRUST REQUIRED / FACTOR C 


With the above equation, it 1s possible to find the required fuel flow for 


different specific conditions without interpolating the flight manual data. 
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Iii. RESULT 


A. BASE LINE CONDITION 


Given condition 


Gross weight 40,000 lbs 
Drag Index 20 
Altitude Standard Sea Level 


1. Converting velocity versus fuel flow to velocity versus drag 
Military thrust (at sea level) 21,800 Ib 
Fuel flow (at military thrust) 25,680 lb/hr 


Thrust = Fuel Flow x Factor C 
C = Military thrust / Fuel Flow 


at military thrust 
C = 21,800/25,680 = 0.84891 


The thrust required for each velocity can be computed by fuel flow times 
factor C(0.84891). The Table 11 shows the result. 


TABLE 11 
THRUST REQUIRED FOR EACH SPEED 


VELOCITY(KTAS) FUEL FLOW(LBS/H) THRUST REQUIRED(LBS) 
360 Les 6,558 
400 9,014 1 0o2 
440 10,623 9,018 
480 12,462 Ory 
20 14,652 12,438 
560 WBE 14,752 
600 Ze so 18,160 
MIL 25,680 21,800 


2. Computing the minimum thrust and maximum endurance speed. 
Maximum endurance speed can be obtained from maximum endurance Mach 
number graph in NAPTOPS FLIGHT MANUAL (see Figure 2.7). The maximum 


endurance Mach number is M 0.378. And maximum endurance fuel flow can be 


DJ). 


obtained from the maximum endurance fuel flow graph in NATOPS MANUAL (see 
Figure 2.8). 


The minimum thrust can be obtained from the maximum endurance fuel flow 
by multiplying with the factor C computed in step I. 


maximum endurance fuel flow ; 5,600 lb/hr 
Tenin = 2 O00rROreS391 = 7522 alb 
3. Computing the K, 


At maximum endurance speed, the thrust required is equal to the total drag, 


that is twice the induced drag or parasite drag, i.e. induced drag and parasite drag are 
the same at this point. 


al 


min ~ Protal = 2Pi = 2D5 
D, = Cry, (1/2)pV*A 


where Dp, = 4,753.9 { 2 = 2,376.9 Ib 


P = 0.0023769 lb sec”/ft4 
V = 0.378 x 1,116.4 ft/sec 
A = 530 ft” 
Cp, = 2Do/ (pVA) 
= (4,753.9)/(0.0023769 x 422.187 x 530) 
= 0,02117183 


D, = (0.02117183)(1/2)(0.0023769)(530)V~ 
= 0.01333 V7 


Thus K, = 0.0133333 


4. Computing the K, 


Use the same concept as for Ky 
Tin = Ptotal = 2D) = 2Do 


Dy = 2W? / (pAV2)(1/meAR) 


1/meAR 


pa 2 


(0.0023769)(530)(422.18)7(2376.9)/(2 x 400002) 
= 0.16677964 


a2 


D; = (2W?/pA)(1/meAR)(1/V7) 
= (2 x 40,0002 / 0.0023769 x 530)(0.16677964)(1/V2) 


Thus 
Ky = 4.2365 x 10° 


Dy ota) = Phrust required = een” =F K»/V- 


These K; and K, values are modified a little to fit the curve of hand book 
data by the trial and error method introduced in Chapter 2. Thus the thrust required 


equation can be written as 
D = 0.01614V2 + 2.0x108/v2 


But if there were no consideration about the Mach effects, the drag at high 
Mach number, Le. greater than M 0.8, will be much less than the hand book data. Ch 
becomes a function of Mach number above Mach numbers of about M 0.8. In this 
region, the parasite drag curve quite closely resembles a sin? curve, and for 


approximation analysis may be represented by the equation 


- 3 
C =C + (C -C sin M - 0.8)/0.8 
Do Doo ( Dom Doo) \( isc 


Thus the drag equation between the Mach number M 0.8 and M 1.2 


CO 
| 


o = W2pSCp y? 


‘a5 2 
D. = 1/2 pS[C =p iC -C n-{(M - 0.8)/0.8}] V 
= S| Doo ( De Daw {( )/0.8}] 
ae} 2 
Dee iKae-t K1/6 C <2 © M - 0.8)/0.8}] V 
ae eS 1/ Doo § Dy Daw sin” {( )/0.8}] 


ees, 2 
=K,fl + (C -C Cc sin-{(M - 0.8)/0.8}] V 
pL ea Dom Da?) Dan {( )/0.8}] 
where C -C = 0.55 
Dom Doo 


00 
D, = Kj[I + 0.55/0.018 sin? {(M - 0.8)/0.8}]V7 


oy: 


5. Converting to fuel flow from drag 
Fuel flow for each velocity can be computed from thrust required, 1. total 
drag, devided by factor C. 


C = MILITARY THRUST / FUEL FLOW (at military thrust) 


Military thrust and fuel flow can be computed with the equation that found in 
Chapter 2. 


Military thrust = 21,800 - 0.4 x Altitude 
= 21,800 -0.4x0 
= 21,800 
Fuel flow; 
Fy, = 26,000 - 0.64 x Altitude 
= 26,000 - 0.64 x 0 
= 26,000 
F. = F, - 16.65 x Drag Index 
= 26,000 - 16.65 x 20 
= 25,667 
F3, = F, - 0.002 x (Gross weight - 40,000) 
26,667 - 0.002 x (40,000 - 40,000) 
25,667 


Thus C = 21,800/25,667 


0.84934 


Fuel flow = Thrust required / 0.84934 


Table 12 and Figure 3.1 show the results of the calculation and deviation from 
hand book data. 


B. VARIATION NUMBER 1 - GROSS WEIGHT 
Given condition 
Gross weight 50,000 lbs 
Drag Index 20 
Altitude Standard Sea Level 
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The gross weight is changed only from the base line condition. It is necessary to 


consider the change of K, and K>. 


Ki = te SD) © 
As you can see from the above equation, K, is independent of gross weight. 


Thus there was no change in K). 
But K is a function of gross weight (as discussed in chapter 3). 


Ky = 2W2/(p S me AR) 


It can be seen that the K-, is proportional to the squared gross weight, but the 


effects of varing gross weight become more largely. 


Thus K5 = K, (50,000/40,000)2-7 
50,0001b 740,000Ib 


= 2.0x108 x (1.25)! 
= 3.653 x 108 


D = 0.0164V2 + 3.653x108/V2 
Military thrust = 21,800 Ib 


Fuel flow; 
Fy = 26,000 - 0.64 x Altitude 


= 26,000 
F. = F, - 16.65 x Drag index 
= 26,000 - 16.65 x 20 
25,667 


F, = Fs - 0.002 x (Gross weight - 40,000) 
= 25,667 - 0.002 x (50,000 - 40,000) 
= 25,647 
C = 21,800/25,647 
= 0.85 


Fuel flow = Thrust required / 0.85 


Table 13 and Figure 3.2 show the results of calculation and deviation from hand 


book data. 
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C. VARIATION NUMBER 2 - DRAG INDEX 
Given condition 
Gross weight 40,000 Ibs 
Drag Index 40 
Altitude Standard Sea Level 


If the drag index was changed from base line condition, what will be effect on the 


result? It is necessary to consider the K, and K> equation. 


K, = 1/2 pSCp 


2W2 / (pSmeAR) 


K 


As can be seen, K> is a function of gross weight, altitude, wing area and aspect 
ratio, not of drag index, thus K, was unchanged. But K, is a function of altitude, wing 
area and drag index. Thus the appropriate K, value can be computed with equation 

“Ipi4o~ “Ipr20* “Ara! “1207 

Now it is necessary to find the relationship of K, between drag index 20 and 
drag index 40. Furthermore, if the relationship of K, between any drag index and 
reference drag index, the various K, for different drag index can be computed. To 
figure out the relationship of K, for each drag index, let us compute each K, value. 


The computing method is the same as introduced in the base line condition step. 


Z 2 
Do = 1/2 pSCp V 
Cp. = 2D ,/(pSV*) 
Ky = 1/2 pSCp 


For example 
when DI = 0 
Do = Dipl 2 AAAS eee 


Ch = 2(2,221.9)/(0.0023769x530x438.97) 
O 
= (7.018132 
Thus 
Ky = 1/2 Doe 
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1/2 (0.0023769 x 530 x 0.018132) 
0.011534 


Table 14 and Figure 3.3 show the result and the relationship of K, for different 


Cy Ky 
FUEL FLOW KTAS MINIMUM O 


Te [ewes sae [ae | eae Jooslooe 
Pe [esses [soe [ae [aa posse 
Pee [ons [ease [aa | sae leona 
Pe [oes | come [ae | ear conf 


The K, of the base line condition was modified to fit the drag index. As you can 


see in Figure 3.3 and the investigation in Chapter 2, the relationship between K, and 


diag Imdex. 


TABLE 14 
KI VALUE VERSUS DRAG INDEX 






MAXIMUM ENDURANCE THRUST (LB) 














drag index is linear and the equation of Ky, function of drag index 1S 


K, = 0.0135 + 0.000132 x Drag Index 
Thus 
= 0.0135 + 0.000132 x 40 
0.01878 


Ky 
D140 


Thrust required = 0.01878V* + 2x108/v2 


Military thrust = 21,800 lb 
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Fuel flow; 
26,000 


Py 
(Fy + 333.3) - 22.1 x Drag Index 


= (26,000 + 333.3) - 22.1 x 40 
= 25,449.3 


£9) 


F3 = F, - 0.002 x (40,000 - 40,000) 
= 25,449.3 


C = 21,800/25,449.3 
= 0.8566 


Fuel flow = Thrust required/0.8566 
Table 15 and Figure 3.4 show the results of calculation and deviation from hand 


book data. 


D. VARIATION NUMBER 3 - ALTITUDE 


Given condition 
Gross weight 40,000 lbs 


Drag Index 20 
Altitude 4,000 ft 


If the altitude will be changed only from the base line condition, how will it affect 


the K, and K,? 
The equations of K, and K, are 


K, = 1/2 pSCp_ 


2W2/(pSmeAR) 


S) 
As can be seen in the above equation, the K, and K, both are functions of the 


air density, and air density is a function of altitude. Thus the K, and Ky must be 


changed with varying altitude. 
The useful relationships came from the ICAO Report written as 


p = P/RT 
T/T, = 1 - 6.875x10°°H 
P/P,, = (1 - 6.875x10°°H)>-2°0! 
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T = 518.688 x (1 - 6.875x10°°H) 


And the relationship between density of air and vertical displacement is 


p/P, = (P/RT)(PQ/RT 9) 
= eae 


Since the effects of changing altitude become larger, it is necessary to modify the 


equation. 


Let Ko = (p/P, 
Then K, is proportional to the Kg and K> is inverse proportional to the Ko. 


yl.2 


Thus in this case 
Dip BR 1G) ©) 


= (1 - 6.875x107°x4,000)>:2°61/(1 - 6.875x10°© x 4,000) 
= 0.88809 


K, = K, . x (0.88809)!-2 
l les 
= 0139976 


2 1.2 
Ky = Ko, / (0.88809) 


2.306x10° 


Thrust required = 0.0139976V2 + 2.306x108/v2 


The military thrust will be reduced as discussed in Chapter 2. 


1 $/0(1/2) 


mes 


aval 
& = Pg ooor / Psst 


1 ,827.69/2,116.22 
0.8636578 


6 = P4 oooft / Pssl 


0.0021 109/0.0023 769 
0.8880895 


= 21,800 x 0.8636578/0.8880895!/2 
19,779 


avail 
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Fuel flow; 
F, = 26,000 - 0.64 x Altitude 
= 26,000 - 0.64 x 4,000 
= 23,440 
Fo = F, - 16.65 x Drag Index 
= 23,440 - 16.65 x 20 
23,107 
F, = F, - 0.002 x (Gross weight - 40,000) 
= 23,107 - 0.002 x (40,000 - 40,000) 
= 0.86463 


C = 19,979/23,107 
= 0.86463 


Fuel flow = Thrust required/0.86463 


Table 16 and Figure 3.5 show the results of calculation and deviation from hand 


book data. 
E. VARIATION NUMBER 4 - WEIGHT, DRAG INDEX AND ALTITUDE 


Given condition 
Gross Weight 50,000 Ibs 
Drag Index 40 
Altitude 8,000 ft 


If ell conditions are changed from base line conditions, it is neceassary to 


consider for each changed condition step by step, same as previous sections. 
At first, think about gross weight change. The changing gross weight affects the 


K» only, because K, is independent of gross weight. 
Ky = 2W? / (pSmeAR) 
As you can see in section 3.2 


K, = K, (50,000/40,000)? 7 
50,0001bs 40,0001bs 
2.0x10° (1.25)*:? 


3.653x10° 
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Next, think about the drag index change. If the drag index is changed, it is 


necessary to consider K, only, because the K4 is independent of the drag index change. 
K, equation, function of drag index is 


K, = 0.0135 + 0.000132 x Drag Index 
Thus | 


Ky)... = 0.0135 + 0.000132 x 40 
D140 


= 0.01878 
And finally, think about the altitude change. 


As you can see in section 3.4, the altitude change affects both K, and K5, 


because K, and K> are function of air density. The equation in section 3.4 can be used 
in this step. 


p/P. = (P/RT)P,/RT,) 


=(P/P,\(T,/T) 
Let Ko = (p/p,)! 

K, = f (p) 

K» = f(1/p) 


: 1.2 
Kanak 
1 = BL yap * §P8,000ft!/ PssP 
= 0.01878 x (0.0018683/0.0023769)!- 
= 0.01406755 


1.2 
K,=K 
= 3,653x108 / (0.0018683/0.0023769)!:2 
= 2.9157x108 


Thrust required = 0.01407V2 + 2.9157x108/V2 


The military thrust will be computed as discussed in Chapter 2. 


7 1/2 
Tg oooft = Fssy (0/6 e) 


& = Pg ooo / Pss1 
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1571.88/2116.22 
0.742777 


© = Pg oooft / Pssi 


= 0.0018683/0.0023769 
SU oun 


a ’ | 
Tg goof = 21,800 x (0.742777/0.786024!/4) 
= 18,264lb 


Fuel flow; 
F, = 26,000 - 0.64 x Altitude 
= 26,000 - 0.64x8,000 
= 20,880 
(Eyemece) - 22lex Drag Index 
= 21,213.3 - 22.1 x 40 
= 20,329.3 
F, = F, - 0.002 x (Gross weight - 40,000) 
= 20,329.3 - 0.002 x (50,000 - 40,000) 
= 20,309.3 


>) 


C = Muilitary thrust / Fuel flow 
18,264/20,309.3 
0.8993 


I 


Fuel flow = Thrust required/0.8993 


Table 17 and Figure 3.6 show the results of calculation and deviation from hand 


book data. 


VARIATION NUMBER 5 - WEIGHT, DRAG INDEX AND ALTITUDE 


We have seen how the fuel flow was changed with various conditions. Not much 


time is needed to find the fuel flow for conditions which can be read directly from hand 
book, but most of the conditions do not correspond to the published conditions which 
can be read directly from the hand book data sheets. Thus to find the fuel flow for 
specific conditions that are not printed in the hand book, one must spend lots of time 


to interpolate each item. The specific fuel flow will now be computed by both methods. 
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Given condition 
Gross weight 48,000 lbs 


Drag Index 27 
Altitude 11,000 ft 
Velocity 455 KTS 


The procedure to compute the K, and K> is the same as in the previous section. 


Consider each different condition. 


1. Gross weight 


K> = K> x (48,000/40,000)*/ 
48,000lbs 740,0001bs 


2.0x108 x (1.2)?-/ 
3.27x108 


K, is independent of changing gross weight. 


2. Drag index 
K, = 0.0135 + 0.000132 x Drag Index 
= 0.0135 + 0.000132 x 27 
= 0.01706 


K- is independent of changing Drag Index. 


3. Altitude 
1.2 
ink 
eae ieee, Oiircoor 1 Peep 


0.01706 x (0.0017007/0.0023769)! 4 


0.01142 


1.2 
re, ape 
D 248 0001b / (P11 0008 / Psst) 


3.27x108/(0.0017708/0.0023769)!-2 


4.8893x10° 


Thrust required = 0.01142V2 + 4.8893x108/V2 
= 0.01142 x (455x1.69) + 4.8893x108/(455x1.69) 
= 7,572 lb 


Compute the factor C: 
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Military thrust = 21,800 - 0.4 x Altitude 
= 21,800 - 0.4x11,000 
= 17,400 Ib 


Fuel flow; 
PY 


20,840 - 0.37 x (Altitude - 8,000) 

= 20,840 - 0.37 x (11,000 - 8,000) 

= 19,730 
Fo) = (i) % 333.3) - 22 x Dragiindex 

= 20,063.3 - 22.1 x 27 

= 19,466.6 
F, - 0.002 x (Gross weight - 40,000) 
= 19,466.6 - 0.002 x (48,000 - 40,000) 
= 19,450.6 


Py 


C = Military thrust / Fuel flow 
= 17,400/19450.6 
= 0.89457 


Thus 

Fuel flow = Thrust required / C 
7,572/0.89457 
8,464 lbs/hr 


Now the fuel flow for given conditions will be computed by interpolating the 
hand book data. Figure 3.7 shows the whole diagram to compute the desired fuel flow 


of given conditions. 


Computations for the steps are: 
1) Fl = (8873 -7988)x(27 - 20)/(40 - 20) + 7988 = 8297.75 
2) F2 = (10324 - 9223)x(27 - 20)/(40 - 20) + 9223 = 9608.35 
3) F3 = (455 - 440)x(9608.35 - 8297.75)/(480 - 440) + 8297.75 = 8789.23 
4) Gl = (7729 - 6962)x(27 - 20)/(40 - 20) + 6962 = 7230.45 
5) G2 = (8840 - 7933)x(27 -20)/(40 -20) + 7933 = 8250.45 
6) G3 = (455 - 440)x(8250.45 - 7230.45)/(480 - 440) + 7230.75 = 7613.25 
7) L4 = (11000 - 8000)x(7613.25 - 8789.23)/(12000 - 8000) + 8789.23 = 7907.24 
8) H1 = (9239 - 8366)x(27 - 20)/(40 - 20) + 8366 = 8671.55 
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DIAGRAM FOR DESIRED FUEL FLOW 
Mobeowns 8,000 FT moO FT Pane ee UO oO UG Lr 


A/S DRAG INDEX A/S DRAG INDEX 


ktas | 20 | 27 | 40 KTAS 
440 Uecker s: Ei S873 | 440 | 6962 Git ee 


455 F3 455 G3 
9223 | F, |10324 7933] G> | 8840 
GROSS WEIGHT 50,000 LBS 


AEPLIUDE. &, | ALTITUDE 8,000 Fr leat 22,000 Fr | Gone jolg ALTITUDE | ALTITUDE 12,000 PP O00 FT 


| DRAG INDEX INDEX DRAG INDEX 


= cee SccE 
eel. fet et 
480 9522 10688 315] 32 | 9216 


G/WT 40,000 LBS } 48,000 LBs | 000 LBS | c/wr 50,000 LBS 50,000) LES 


DRAG INDEX | DRAG INDEX INDEX 
ZZ 





Figure 3.7 Diagram for Desired Fuel Flow. 


8] 


9) 
10) 
11) 
12) 
13) 
14) 


15) 


H? = (10688 - 9522)x(27 - 20)/(40 -20) + 9522 = 9930.1 
H3 = (455 - 440)x(9930.1 - 8671.55)/(480 -440) + 8671.55 = 9143.51 
J1 = (8192 - 7447)x(27 - 20)/(40 - 20) + 7447 = 7707.75 

J2 = (9216 - 8315)x(27 - 20)/(40 - 20) + 8315 = 8630.35 

J3 = (455 - 440)x(8630.35 - 7707.75)/(480 - 440) + 7707.75 = 8053.73 


M4 = (11000 - 8000)x(8053.73 - 9143.51)/(12000 - 8000) + 9143.51 
8326.174 


De ass (48000 - 40000)x(8326.17 - 7907.24)/(S0000- 40000) + 7907.24 


Desired fuel flow for given condition is 8,242 lbs/hr. 


Thus deviation of the results 


(8,464 - 8,242)/8,242 x 100 = 2.7% 
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lV. CONCLUSTON 


It 1s concluded from the investigations that computing the fuel flow rate from the 
computer program results in less deviations and is simpler than interpolating from 
NATOPS FLIGHT MANUAL DATA. The total lines used in this program is 53 lines 
for the min frame program and 229 lines for the hand-held calculator program, and 
both programs are inserted in Appendix. 

Though the deviation of fuel flow for this program was restricted to a maximum 
of 10 percent, for most conditions deviations are distributed below 5 percent. To 
determine the fuel flow at specific conditions which can not be read directly from the 
flight manual, it will be necessary to calculate fifteen steps of interpolations. This 1s 
very time consuming and the probability of making a mustake is very high when 
calculating so many interpolations. 

Thus the computer program to compute the fuel flow rate can not only save time 
for many conditions which must be changed continuously during flight, but it also 
reduces the probability of making a mistake in interpolations steps. 

It is recommended that a similar computer program be developed for specific 
flight stages, i.e., take-off, climb, or landing, using the same concepts that were used in 


this program. 
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00010 
00020 
00030 
00040 
00050 
00060 
00070 
00080 
00090 
00100 
00110 
00130 
00160 
00170 
00180 
00190 
00200 
00210 
00220 
00230 
00240 
00250 
00260 
00270 
00280 
00290 
00300 
00310 
00320 
00330 
00340 
00350 
00360 
00370 
00380 
00390 
00400 
00410 
00420 
00430 
00440 
00450 
00460 
00470 
00480 
00490 
00500 
00510 
00520 
00530 


APPENDIX A 
COMPUTER PROGRAM (BASIC) 


K2=2.0E+08 ! REM **% BASED ON S.S.L. AND GROSS WEIGHT GOOOOLBS xx 
PRINT ‘INPUT ALTITUDE,GROSS WEIGHT,DRAG INDEX,AIR SPEED,A/S INC’ 
INPUT A,W,19,V1,VP 

K1=(0.0135+.000132*%1I9) ! 
D4=1-6.875c-6XA ! 


REM *® Kl VARIES WITH DRAG INDEX Xx 

REM *% ICAO SCALE BASE Xx 

T1=518.688%D4 ! REM *% ABSOLUTE TEMP (CR) BASED ON ALTITUDE Xx 
D5=D4**5.2561 ! REM XX PRESSURE RATIO BASED ON ALTITUDE xXx 
A9=(1.4%1714.87%T1) !' REM ** SPEED OF SOUND (SQUARED) xx 
K9=(D5/DG)*x*1.2 ! REM ** ALTITUDE CORRECTION FOR 'Q* xx 

K1=K1*¥K9 ! REM ¥% PAXMX¥2 ALTITUDE CORRECTION xx 

K2=K2/K9 ! REM ¥* PAXM**2 ALTITUDE CORRECTION xXx 
K2=K2*XC(W740000)%*2.7 'TREM*X® WEIGHT OTHER THAN GOOOOLB CORRECTIONXX 
PRINT 'K1l=';K1, 'K2=';K2 

PRINT TABC1); 'VELCKTS)';TABC12)3 "MACH; TABC20); "PAR DRAG’; 

PRINT TABC34);*PRO DRAG; TABC48); "MACH DRAG’; TAB(C62);*FUEL FLOW’ 
V2=V1¥1.68894 

M=V27A9%x0.5 


D8=1 ! REM ** CDO MACH CORRECTION BASE Xx 

IF M>.8 THEN 480 

DI=K1¥V2"X2 !' REM ®¥ PARASITE DRAG (CSUBSONIC) xx 
D3=D1XD8 ! REM ®¥ PARASITE DRAG (MACH CDQ) xXx 
D2=K2/V2xx2 !' REM XX INDUCED DRAG Xx 

DO=D3+D2 ! REM *% MACH DRAG Xx 

D7=D1+D2 ! REM ®¥® SUBSONIC TOTAL DRAG xx 


T=21800-AX0.4 ! REM *® MIL THRUST OTHER THAN S.S.L. CORRECTION Xx 

IF A<=8000 THEN 320 

eae bath REM ** F/FLOW CORRECTTION ABOVE 8000 FTX 
> 

F1=26000-0.64000%A ! REM *®X F/FLOW CORRECTION BELOW 8000 FT Xx 

IF I19>20 THEN 360 

F2=F1-16.65*%1I9 ! REM **® F/FLOW CORRECTION BELOW DRAG INDEX 20 x 

GO TO 370 

F2=(F1+333.3)-22.10%19 ! ® F/FLOW CORRECTION ABOVE DRAG INDEX 20 xX 

IF W<=GO0000THEN 400 

F3=F2-C(W-40000)%0.002 ! REM *® F/FLOW CORRECTION WITH WEIGHT xx 

GO TO 410 

F3=F2 

C=T/F3 ! REM ®¥% COEFFICIENT FOR FUEL FLOW Xx 

FF=D0/C ! REM Xx FUEL FLOW Xx 

PRINT TABC3)3;V1;TABC8) 5M; TABC18);D2;TABC32);D3; TABC46); DO; 

PRINT TABC60);FF 

V1l=V1+VP 

IF V1>610 THEN 510 

GOTO 190 

D9=.55XCSINCOM~.8)7.8)**3)7.0180 

D8=1+D9 

GOTO 230 

PRINT "C=";C,*T=';T 

STOP 

END 
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APPENDIX B 


COMPUTER PROGRAM (HAND-HELD CALCULATOR) 


PieLBL *FE 
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a7 STa 24 
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pe -y2?° 
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11 ST) 8 
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‘5 STD #2 
LeeLBL ¥¥" 
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19 PROMPT 

2 37) @3 
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29 ,@a0132 
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a} + Zoi RCL 23 
ioc 1S on? “FFs° 
[oot 23 ARGhos 
154 570 24 26s A lEw 
aCe chilies Sa 2e5 STOP 
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